Boltzmann transport and residual conductivity in bilayer graphene 
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A Drude-Boltzmann theory is used to calculate the transport properties of bilayer graphene. We 
find that for typical carrier densities accessible in graphene experiments, the dominant scattering 
mechanism is overscreened Coulomb impurities that behave like short-range scatterers. We antic- 
ipate that the conductivity (j(n) is linear in n at high density and has a plateau at low density 
corresponding to a residual density of n* — ^yn■nnp^n, where n is a constant which we estimate using 
a self-consistent Thomas-Fermi screening approximation to be n « 0.01 q^p ~ 140 x 10^" cm~^. 
Analytic results are derived for the conductivity as a function of the charged impurity density. We 
also comment on the temperature dependence of the bilayer conductivity. 

PACS numbers: 81.05.Uw; 72.10.-d, 73.40.-c 
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I. INTRODUCTION 

The recent experimental realization of a single layer 
of carbon atoms arranged in a honeycomb lattice has 
prompted much excitement in both the theoretical and 
experimental physics communities (For a recent review, 
see Ref. [l] and references therein). The focus of the 
current work is on bilayer graphene which has received 
less attention both theoretically and experimentally, but 
is nonetheless of equal importance both for technolog- 
ical application and for fundamental science. Bilayer 
graphene is two monoatomic layers of graphene sepa- 
rated by about 0.3 nm, which is the interplane distance in 
bulk graphite. Similar to single layer graphene, bilayer 
graphene has been realized experimentally through the 
mechanical exfoliation of graphite onto Si02 substratesii 
While the band structure of a single layer of graphene 
has a linear dispersion, theoretically bilayer graphene has 
a quadratic dispersion with an effective mass of about 
0.03 rUe making it similar to the regular two dimensional 
electron gas (2DEG). Despite the quadratic spectrum, bi- 
layer graphene shares two important features with single 
layer graphene (hereafter referred to simply as graphene) 
that distinguish it from regular 2DEGs. First, the bilayer 
effective Hamiltonian?.^!^!^!^ is chiral which gives rise to 
the anomalous integer quantum hall effect fi Second, un- 
biased bilayer graphene is a semimetal implying that one 
continuously moves from electron-like carriers for posi- 
tive gate voltages to hole-like carriers for negative gate 
voltages without any gap in the spectrum. We note that 
although recent experiments^ on graphene bilayers have 
been able to open a gap by connecting the upper layer 
to an external top gate, here we ignore this additional 
degree of iieedom^^ 

By considering the gapless bilayer situation, the low 
density transport resembles that of single layer where 
Coulomb impurities in the substrate create an inhomo- 
geneous density profile breaking the system into puddles 
of electrons and holes. The bulk residual density n* in- 
duced by these impurities has been calculated for single 
layer graphen o^^i^^ which shows agreement with recent 



experimental studies J^ii^ii^ The high density transport 
in single layer graphene with screened Coulomb impuri- 
ties was discussed in Refs. The goal 
of the present work is to generalize these high-density 
and low-density single-layer graphene Boltzmann trans- 
port theories to the case of graphene bilayers. We note 
that the Boltzmann transport theory developed here ig- 
nores the effects of phase-coherence which was studied in 
Refs. [mi. 



II. BILAYER HAMILTONIAN AND 
BOLTZMANN TRANSPORT 

The effective Hamiltonian for bilayer graphene is now 
well established in the theoretical literature (See Refs. 0, 
[114,5,6,23,24,25). First principles and band structure cal- 
culations show that at both very small energies and very 
large energies, bilayer graphene has a linear spectrum. 
For energies 2 x 10~^ eV e < 0.1 eV bilayer graphene 
has a quadratic spectrum (see e.g. Refs. y|j). In princi- 
ple, bilayer graphene should have a complicated crossover 
Hamiltonian moving from linear to quadratic and back 
to linear as one increased the carrier density. To our 
knowledge, a theoretical framework to understand the 
structure and effects of chirality within this crossover 
has not yet been developed. However, as we argue be- 
low, charged impurities in the substrate induce a residual 
density n* in bilayer graphene that corresponds to a typ- 
ical Fermi energy e > 0.01 eV which is larger the the 
lower energy scale for using the quadratic Hamiltonian, 
and the range of experimental gate voltages Vg < 50 V 
induces a maximum carrier density with Fermi energy 
e < 0.1 eV which is comparable to the limit where 
the low-energy effective quadratic Hamiltonian begins to 
break down. Therefore, for realistic samples, the condi- 
tion 2 X 10^'' eV < e < 0.1 eV is mostly satisfied and 
the quadratic Hamiltonian Ti. proposed by McCann and 
Falko should be an excellent approximation, where in ad- 
dition bilayer graphene is weakly interacting in this en- 
ergy window. In what follows, we use^ 
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FIG. 1: (Color online) Scattering cross-section as a function 
of angle. Unlike single layer graphene, both bilayer graphene 
and 2DEG are dominated by backscattering. 



TABLE I: Summary of Boltzmann transport results in 2 
d electron gas (2DEG), single layer graphene and bilayer 
graphene. For screened Coulomb scattering results in 2DEG 
and bilayer graphene we assume that gTF/2fcF > 1 (see 
text), and observe that arising from different physics, screened 
Coulomb scattering gives cr ~ n in all three cases. 
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This Hamiltonian can be diagonalized giving e/c = 
±h'^k'^/2m where m = 271^2/(37^0^) « 0.033 rrie, and 
7o ~ 3.16 eV is the in-plane coupling and 71 « 0.39 eV 
is the out of plane coupling, a « 0.246 nm is the lattice 
constant and rrie is the electron mass. The eigenvectors 
^± ~ (e^'2^'',±l) give the aforementioned chiral proper- 
ties where k — keKp(i9k)- 

Using this diagonal basis, one can calculate the scat- 
tering time T using the Boltzmann transport theory^^ to 
find 
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' i-i V /I V ^ K.K. that tnr iinsrrppnpn ( ,mi nmh imnnritips 



Before calculating the Boltzmann conductivity a — 
{2e^ /h)kFVFT, we first use dimensional arguments to de- 
termine the dependence of conductivity on carrier den- 
sity. In 2DEGs and bilayer graphene the Fermi velocity 
vf — hkF/rn depends on carrier density through n ^ k^ 
while the inverse screening length gxF = 4,me'^ / (kIt?) is 
density independent. Since the gas density parameter 
Ts ^ scales as n~^^'^. This is all in sharp con- 

trast to single layer graphene where vf ~ 10^ m/s is 
constant and qtf = 4fcFf's depends on carrier density, 
and rg = /{kvf) ~ 0.8 is a density independent con- 
stant that depends mostly on the dielectric constant of 
the substrate. 

and 



' fc| giv- 



where the momentum transfer q = |k — k'| = 
2fcF sin(0/2), v{x) is the screened scattering impurity po- 
tential and k • k' = cos 61. The effects of chirality are 
captured in F^^^,. For a regular 2DEG whose Hamilto- 
nian is not chiral, we have F{6) = 1, whereas F{6) ~ 
(1 -f cos6')/2 for graphene and F{e) = (1 + cos26l)/2 for 
bilayer graphene. The term (1 — k-k')Fi^^, determines the 
scattering cross-section which has been plotted in Fig. [21 
We observe that in contrast to single layer graphene, 
both bilayer graphene and the 2DEG are dominated by 
backscattering {6 = tt). Introducing x — q/(2kF), we 
find 



l^J^dx\v{x)\ 
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for bilayers, 



"•imp||T dx\v{x)\^x'^\/l - x^ for graphene. 



(3) 



where the single la yer gr aphene result was reported pre- 
viously in Refs. [T]]fl9ll27l . 



that for unscreened Coulomb impurities tq 
ing (Tc ~ n^, whereas overscreened Coulomb scatterers 
behave similar to white-noise disorder giving density in- 
dependent T and a ^ n/riimp (this is similar to a low- 
density 2DEG, where Coulomb scatterers are strongly 
screened if qtf ^ 2A;f). Writing (j(n) ~ n", we note 
that a = 1 in both the linear and quadratic Hamiltonians 
arising from very different reasons (See Table The so- 
lution within a crossover between a quadratic and linear 
Hamiltonian (see above) is beyond the scope of this work, 
but we observe that approaching the crossover from ei- 
ther side gives a > 1, and for Coulomb scatterers located 
at the Si02 interface, we have a < 2. In what follows we 
focus on the experimentally relevant regime, where we as- 
sume that qtf /2kF > 1 which is typically called the low 
density regime in 2DEG literature^SS, In this context, even 
in GaAs heterostructures (where m « 0.07 rrie) moving 
to higher density results in a complicated crossover where 
the exponent a slowly decreases with increasing density 
as other scattering mechanisms come into play.^* 

In single layer graphene, it was shown by Refs. H^HO 
that for q < 2/cf the static dielectric function calcu- 
lated in the Random Phase Approximation is identical 
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to that of the much simpler Thomas-Fermi approxima- 
tion. A similar result holds for 2DEGs. While it is not 
clear if this holds in bilayer where the polarizability has 
only been calculated numerically (See Ref. 31), these re- 
sults indicate that Thomas-Fermi approximation (which 
allows for analytical results) should capture both qual- 
itatively and quantitatively the transport properties of 
bilayer graphene. Within the Thomas-Fermi approxima- 
tion, the potential of a charged impurity located at a 
distance d from the substrate is 



v{q) 



~qd 



2m 



(4) 



where in the second equation we have used the further 
approximation (also called "complete screening approx- 
imation") that qtf = ^me^ / [nh?) ^ 1 nm"-"^ is larger 
than the maximum transfered momentum q < 0.3 nm^^. 
Herein lies an important difference between single layer 
graphene and bilayer graphene. For single layer graphene 
gxF — 4fcF7's depends on density, so that both the 
screened and unscreened Coulomb potential scale as k^^ . 
It is this property of single layer graphene that gives rise 
to the conductivity with Coulomb scatterers being linear 
in density and the inapplicability of Gaussian white-noise 
models (i.e. zero-range scattering) to capture the trans- 
port properties. In contrast, for bilayer graphene and 
2DEG, gxF is a density independent constant which is 
larger than the typical momentum transfered in current 
experiments, and therefore the strong screening approx- 
imation ((7tf > 2k-p) applies except at very high car- 
rier densities. In this context, bilayer graphene is much 
more similar to 2D Si MOSFETs than to single layer 
graphene. Analytic results for the conductivity can be 
obtained both in the limit d — > and 2A;f/9tf ^ 0. 
Keeping both to leading order, we find 
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(5) 



The result a = {Ae^ / ■nh){n / rii-^-p) , where conductiv- 
ity is linear in density is valid for d~^,q'j:Y S> 2fcF, al- 
though we expect deviations from this linear behavior for 
n > 6 X 10^^ cm~^, a regime, which in principle should be 
accessible in future experiments. We note that the lin- 
ear in density behavior was anticipated in Ref. and in 
Ref. [24b ut we point out that the low density saturation 
in Ref. 5' arises from a completely different and universal 
mechanism^"^'^^ that we believe is unobservable in current 
bilayer graphene samples because of the large and non- 
universal n* arising from the disorder induced electron- 
hole puddles liii^ In the following section we calculate 
the voltage fluctuations and residual density induced by 
charged impurities. 



III. LOW DENSITY RESIDUAL DENSITY 



This section follows closely the derivation in Refs. [Tlll3^ 
and applying the same formalism to the case of bilayer 



graphene. We consider iVimp impurities located at the 
points {n} in a 2D plane. If each impurity has a potential 
4>{r), then the disorder averaged potential 

V ^ dr,dr2 ■ ■ ■ P[ri]P[r2] • • • ^ cj){n), 

= nimp4'{q = 0), (6) 

where to get the second line we have assumed that the im- 
purities are uncorrelated and uniformly distributed. For 
sample area A, rii^p = Ni^p/A, P[ri] = A~^ and 4'{q) is 
the 2D Fourier Transform of the impurity potential 0(r). 
For example, the real space Coulomb potential for an im- 
purity located a distance d from the graphene (or bilayer) 
plane, (j}{r) = {e^ / k)[P + d"^]^^^^. This gives for the 
bare potential 4)o{q) — {2TTe'^ / k) exp{—qd)/q. Using the 
Thomas- Fermi dielectric function e(g) = 1-1- qtf/? gives 
the screened Thomas-Fermi potential shown in Eq. |4l 
Since the static polarizability at g ^ is related to the 
density of states v by the compressibility sum rule, the 
result for the disorder averaged potential V = rii^p/v 
is actually quite general. For bilayer graphene v is con- 
stant, and the threshold voltage shift can be immediately 
obtained from n — Vv — riimp- It is the property that 
qxF is independent of carrier density that explains why 
the threshold voltage shift in bilayers is n = nimp, which 
is in contrast to graphene^i where the density dependent 
inverse screening length gives n — n?^-^p/4n*. This non- 
linear dependence of the threshold voltage on charged im- 
purity density in monolayer graphene has recently been 
verified experimentallyj^^ 

The disorder averaged potential fluctuations (U^) are 
obtained in a similar fashion. 



(O = / dridr2---P[ri]P[r2]' 



■E 
■E 

1=3 = 1 



dq 

(2^ 



(7) 



For the Thomas- Fermi potential Eq. [4l we have 



(U2> - U2 



2tt 
— 27m 



qdq 



27re2 
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-I 2 



K q 



imp 



9tf_ 
= 2gTFd), 



(8) 



where Cq^(x) — dx{xe^ Ei[x\) and Ei[x\ = dte^*- /t is 
the exponential integral function. While the results pre- 
sented here are for the Thomas-Fermi potential Eq. IH 
it is straightforward to generalize this result to obtain 
(^RPA Jqj. jipa screening (e.g. using the numerical bi- 
layer dielectric function calculated in Ref. Isil) . but for the 
relevant density scale set by (j'TFj we expect these results 
to be quantitatively quite similar. The self-consistent 
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FIG. 2: (Color online) Self-consistent Boltzmann theory for 
bilayer graphene (solid lines) compared with results of Ref. [ill 
for monolayer graphene (dashed lines). 



density n* is found by setting Ep — {SV"^) (see dis- 
cussion m Refs.[ilil) where this approximation ignores 
the exchange and correlation contributions which are be- 
Heved to be smalL— i^Sii^ Applying this self-consistent 
procedure to bilayer graphene, we find (nn* /2m)'^ — 



and 



q^pCr{2qTFd) 



2-K 



(9) 



140 X 10^°cm"2 



where we approximate C'^ {x 



1) 



0.085. This disor- 



der induced residual density gives rise to a non- vanishing 
conductivity even as the external gate voltage is tuned 
through zero. 



IV. SELF-CONSISTENT TRANSPORT THEORY 

Combining the low and high density results developed 
in the preceding sections, we find for bilayer graphene 



cr(n 



n) 




n < n 



n > n 



(10) 



where n — 



and 



140 X 10i*'cm-2. This re- 



sult predicts that a reasonably clean bilayer sample with 
JT-imp = 5 X lO^^cm"^ would have a mobility /j, ^ « 
6000 CTo^/Vs. The residual density n* = ^ ^imp" ^ 
25 X 10^°cm^^ with a plateau width ^ AV and mini- 
mum conductivity cr^m ~ 7e^//i. When compared to re- 
cent experimental results)^ these estimates agree well for 
the mobility, plateau width and minimum conductivity, 
although not for the offset gate voltage determined from 
fi (see Refs. [ssMsgI for a discussion of other factors that 
could determine the threshold voltage shifts and could 



account for this discrepancy). These results do not de- 
pend qualitatively on the precise choice of d, although 
the results do depend quantitatively; for example, for 
the same value of Tiimp, increasing c? by a factor of 2 gives 
n* « 15 X 10^°cm^^. This may be important for bi- 
layer graphene, since the distance between the two layers 
c ~ 0.3 nm suggests that for the same substrate, the effec- 
tive distance from the charged impurities would be larger 
for bilayers than for graphene. The results of Eg . [TOl are 
shown in Fig. [5] for both a clean (njmp = lO^^cm"^) and 
dirty (nimp = 5 x lO-'^^cm^^) samples and compared with 
the results of Ref. [ll| for graphene using the same charged 
impurity densities and keeping d = 1 nm fixed. One 
notices immediately that for the same charged impurity 
concentration, graphene has a factor of 16 higher mobil- 
ity, smaller plateau widths and larger minimum conduc- 
tivities than the bilayer system. These predictions can 
be easily tested in future experiments. We also note that 
only a(n) ~ n behavior was observed in the experiments 
of Ref. |3J indicating that the current experiments are 
adequately described by the complete screening approx- 
imation, although future experiments should observe a 
"super-linear" conductivity (i.e. a > 1) at higher den- 
sity arising both from the high-density Thomas-Fermi 
corrections as well as from the high-density graphene bi- 
layer Hamiltonian crossing over from quadratic to linear. 



V. TEMPERATURE DEPENDENCE 

Shown in the Fig[l]is the effect of chirality on the domi- 
nant scattering angle, where the suppression of backscat- 
tering seen in graphene is absent for bilayer graphene. It 
was argued recently that the suppression of 2fcF scatter- 
ing in graphene implied weak temperature dependence^ 
until higher temperatures where phonon effects are ob- 
served. The fact that 2kp scattering is not suppressed 
in bilayer graphene does not, however, necessarily lead 
to any screening (or equivalently, Friedel oscillation) in- 
duced strong temperature dependence in the resistiv- 
ity. The temperature dependence in bilayer graphene 
depends on three dimensionless parameters: gTr/Sfep; 
T/Tf; and T/Tu where Td « h/2T is the Dingle tem- 
perature and Tp is the Fermi temperature. The temper- 
ature dependence from screening will be weak if any one 
of these three parameters is not large. 

The actual value of the dimensionless screening param- 
eter 



gxF 2me 
2me^ 



(11) 



is reasonably small even at the lowest accessible carrier 
density set by n* ~ 2.5 x 10"'^"'^ cm~^, making the tempera- 
ture dependence arising from screening rather weak. Sec- 
ond, the dimensionless temperature T/Tp is rather small 
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since the Fermi temperature Tp changes from 120 K 
at low density (set by n*) to 1200 K at high density 



(set by n = kp/ir = aVg 



3.6 X 10^^ cm ^, where 



■ - ■ T / y — 7 ■ ■ 

7.2 X lO^*' cm^^V^ is a geometry related factor). 

The Fermi temperature is relatively high due to the 
very small carrier effective mass (« 0.03 rrie) in bi- 
layer graphene. In fact, the effective mass for bilayer 
graphene is less than half of that for 2D electrons in GaAs 
(~ 0.07 TOe), where the temperature dependence arising 
from screening is extremely small^^ even at much lower 
carrier densities. Therefore, we do not anticipate any 
strong screening-induced temperature dependence in bi- 
layer graphene resistivity. Finally, the currently available 
bilayer graphene samples have very small mobilities re- 
sulting in relatively strong collisional broadening effects 
(i.e. high Dingle temperature) which would suppress any 
small screening induced temperature dependence that 
could have arisen at low temperatures. In particular, 
a mobility of 5000 crn^/Vs as observed in Ref. [13 corre- 
sponds to a Td ~ 50 if, leading to further suppression 
of any screening induced temperature dependence in the 
conductivity. 

It is therefore gratifying to see that the recent ex- 
periment on bilayer graphene^i does not see much tem- 
perature dependence in the low temperature resistiv- 
ity in spite of the importance of 2fcF scattering in bi- 
layer graphene. The temperature dependence seen in the 
plateau region is likely to be caused by thermal popula- 
tion of carriers since Tp ~ 120 K, and the temperature 
dependence is seen for T > 100 K. For T T^, the 
thermally excited carrier density n{T) = J u{e)f{e, T)de, 
where f{e,T) is the Fermi distribution function. This 
gives n ^ T for bilayer graphene, while n ^ 
for graphene giving for the conductivity (ignoring any 
phonon or electron-hole scattering contributions) 



8 In 2 me'' 

— W^(fcr)2 



-{kT) for bilayers. 



for graphene. 



(12) 



Note that the thermal excitation of carriers leads to 
an enhanced^ cr(2^) whereas temperature dependent 
2kp screening typically suppresses cr{T). For bilayer 
graphene, this result suggests that for T « 260 K, there 
would be a 300 percent enhancement in the minimum 
conductivity which is consistent with the observations 
and estimates of Ref. 34. The situation for single layer 
graphene is quite different, since even in the low density 
saturation regime, T/T-p <C 1. For example, a residual 
density n* — 2.5 x lO^^cm^^ corresponds to Tp ~ 100 K 
for bilayer graphene and Tp « 700 K for single layer 
graphene. For this reason, one expects the minimum con- 
ductivity in bilayer graphene to show stronger tempera- 
ture dependence than single layer graphene even though 
Eq. [13 shows the bilayer conductivity scaling as ^ T 
(compared to cr ~ for the monolayer) . 



VI. CONCLUSION 

The formalism developed here captures many features 
of graphene bilayer transport. In particular, we show 
that both the low density saturation and the high-density 
linear in density behavior seen in recent experiments 
arise from the same charged impurities that are invari- 
ably present in samples exfoliated onto a Si02 substrate. 
We note that the agreement between our self-consistent 
Boltzmann theory and available graphene (and bilayer) 
transport experiments and direct measurement of the 
electron and hole puddles giving rise to the residual den- 
sity, indicate that in current experiments, most of the 
charged impurities reside close to the graphene-substrate 
interface (i.e. within ^ 1 nm) similar to the correspond- 
ing situation in Si-Si02 MOSFET structures. If the typ- 
ical distance of impurities from the graphene (or bilayer) 
sheet could be changed in future experiments (e.g. by 
removing the Si02 substrate), the formalism we have 
developed predicts quantitatively its effect on transport 
properties and on the magnitude of the electron and hole 
puddles. Another conclusion of our work is that bilayer 
mobilities are an order of magnitude smaller than cor- 
responding mobilities for graphene samples on similar 
substrates. We also argue that the minimum conduc- 
tivity, plateau width and threshold voltage have differ- 
ent scaling in graphene and bilayers and have different 
dependence on the substrate dielectric constant. These 
differences could potentially be useful when designing a 
particular device application. Moreover, a systematic 
study of the differences in transport properties between 
single layer graphene and bilayer graphene could verify 
our claim that the underlying mechanism for most of the 
observed transport properties on current graphene (and 
bilayer) samples is dominated by charged impurities. 

In addition we have argued that both graphene and bi- 
layers have weak temperature dependence making these 
materials have among the highest room temperature mo- 
bility for any field effect device. One important conse- 
quence for technology is that since mobility is limited by 
charged impurities, better samples can be made either 
by removing charged impurities in the Si02 substrate or 
by using different substrate (e.g. vacuum for the case of 
suspended graphene and bilayers). Attaining high mo- 
bility is also necessary to access new physics such as the 
fractional quantum Hall effect or the "universal" mini- 
mum conductivity, both of which we believe is not seen in 
current experiments due to the large number of charged 
impurities. 

In one sense, bilayer graphene can be thought of as a 
new material that shares some properties with regular 
2DEGs (e.g. quadratic Hamiltonian) and some proper- 
ties with graphene (e.g. chiral Hamiltonian). However, 
theoretically, bilayer graphene is a far more interesting, 
where as discussed earlier, modestly increasing the back 
gate voltage from that of current experiments should in- 
duce a sufficiently large carrier density to see several 
interesting effects. For example, a "super-linear" con- 
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ductivity is predicted both by high-density corrections in 
the Thomas-Fermi approximation (within the quadratic 
Hamiltonian) as weU as by the theoretically expected 
crossover between the quadratic and linear Hamiltoni- 
ans. At even higher (and perhaps realistic) densities, one 
would expect a strong increase conductivity as one pop- 
ulated higher bands. The theoretical framework for un- 
derstanding these high-density effects (including the role 
of chirality) and their observation remains an exciting 
avenue for future theoretical and experimental research. 
In summary, we have proposed a simple theory for bi- 



layer graphene transport including the effects of screened 
Coulomb impurities. The result of our self-consistent 
Drude-Boltzmann semi-classical diffusive transport the- 
or3*ii is in good agreement with the recent experiments 
of Ref. [13 and we make several predictions that can be 
tested in future experiments. 

We thank Andre Geim for sharing with us his unpub- 
lished data (Ref. Is^ ) and for a careful reading of our 
manuscript. This work is supported by U.S. ONR. 
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